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Abstract. We analyze coherently the stellar and nebular 
energy distributions of starbursts and Hii galaxies, us- 
ing our evolutionary synthesis model, PEGASE (Fioc & 
Rocca-Volmerange 1997, 2000), coupled to the photoion- 
ization code CLOUDY (Ferland 1996). The originality of 
this study is to relate the evolution and the metallicity 
of the starburst to the past star formation history of the 
host galaxy. Extinction and geometrical effects on emis- 
sion lines and continua are computed in coherency with 
metallicity. We compare our model predictions to an ob- 
served sample of « 750 H ii regions and starbursts. 

When fitting [O iii]a4363/[0 iii]a5007, [O i]A630o/Ha, 
[Sii]AA67i6,673i/Ha, [N ii] A6584/Ha and [Oiii]a5007/H/3, 
the most striking feature is the decreasing spread in U 
with increasing metallicity Z. High-[/ objects systemat- 
ically have a low metallicity while low levels of excita- 
tion happen at any Z. The best fits of emission line ra- 
tios are obtained with a combination of a high- and a 
low-ionization components. No additional source of ion- 
izing photons - shocks or hidden AGN - is needed. The 
high level of excitation observed in metal-poor H ii galax- 
ies requires a very young population (< 3Myr), while 
starburst nuclear galaxies (SBNGs) are consistent with 
a wider range of age (< 5Myr). 

Colors {B — V, V — R) and equivalent widths are fitted 
in coherency with emission line ratios. An underlying pop- 
ulation is needed, even for small-aperture observations. 
This evolved population not only reddens the continuum 
and dilutes the equivalent width of the emission lines, but 
also participates in the ionization process. Its main ef- 
fect on line ratios is to maintain a high level of excitation 
when the burst stops. Models combining underlying pop- 
ulations typical of Hubble sequence galaxies and instan- 
taneous starbursts with ages between and 8Myr agree 
satisfactorily with all the data. 

Key words: Galaxies: starburst - Galaxies: evolution - 
Galaxies: stellar content - ISM: Hll regions 
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1. Introduction 

The "starburst" phenomenon calls to mind a class of ob- 
jects dominated by the radiation of massive stars (Gallego 
et al. 1995) embedded in a dusty Hii region. ModeUing 
starbursts consistently is the key to interpret the prop- 
erties of the actively star-forming galaxies detected at a 
redshift z > 2 (Giavalisco et al. 1996; Madau et al. 1996; 
Steidel et al. 1996; Lowenthal et al. 1997). To this purpose, 
the extensive study of local samples is a requisite step. 

Local starbursts span a wide range of types, from in- 
dividual H II regions in spiral arms to blue compact H ii 
galaxies and huge kpc-scale nuclear starbursts (Coziol et 
al. 1998). Types differ from one another in their main 
characteristics (line ratios, equivalent widths, colors), and 
presumably have different stellar populations and physical 
conditions. So, is it possible to find correlations between 
the basic properties: age, initial mass function, metallicity, 
relative distributions of stars, gas and dust? Among these, 
which one dominates spectral properties? In which range 
do parameters vary? Answering these questions requires 
a consistent model of starbursts coupling the evolution of 
stars, dust and gas. 

The first studies dealing with emission lines of H ii re- 
gions (Shields 1974, 1978; Stasihska 1978, 1980; McCaU et 
al. 1985) relied on single-star photoionization models. Av- 
erage physical properties of large starburst samples were 
deduced from these pioneering works. As an example, a re- 
lation between the metallicity Z and the ionization param- 
eter U was derived by Dopita & Evans (1986) by fitting 
the emission line ratios from large samples of extragalactic 
Hii regions. 

More realistic stellar populations of star clusters were 
computed by McGaugh (1991) with the Salpeter (1955) 
initial mass function (IMF), to calibrate metallicity de- 
pendant line ratios. Further improvements were the imple- 
mentation of theoretical tracks of massive stars, allowing 
to follow the evolution in time of a starburst (e.g. Garcia- 
Vargas & Diaz 1994). Thanks to the computation of evolu- 
tionary tracks for various metal abundances, the influence 
of the metallicity could also be studied (Cid-Fernandes 
et al. 1992; Cerviiio & Mas-Hesse 1994; Olofsson 1997; 
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References 


Starburst type 


Balmer absorption 


Reddening 






correction 


correction 


Terlovich ot al. 1991 


HiiG, SBNG 


no 


no 


Veilleux & Osterbrock 1987 


SBNG 


no 


no 


Veilleux et al. 1995 


IRAS SBNG 


yes 


no 


French 1980 


HiiG 


no 


no 


Leech et al. 1989 


IRAS SBG 


yes 


yes 


Storchi-Bergmann et al. 1995 


HiiG, SBNG 


no 


no 


Contini et al. 1998 


IRAS SBNG, Hii 


yes 


no 



Table 1. Data sources. SBNG: starburst nucleus galaxies; HiiG: Hii galaxies; Hii: 
H II regions. 



Garcia- Vargas et al. 1995a, b; Stasinska & Leithercr 1996). 
Evans (1991) emphasized the importance of using up-to- 
date stellar atmosphere models to calibrate nebular diag- 
nostics. The most recent models take into account mod- 
ern, updated stellar physics (Leitherer & Heckman 1995; 
Leitherer ct al. 1999), in particular the impact on stel- 
lar spectra of line blanketing and of departures from local 
thermodynamic equilibrium (Stasinska & Schaerer 1997). 

Till now, only a few large datasets covering a wide 
range of physical conditions have been analyzed with 
state-of-the-art models, coupling evolving stellar popula- 
tions and photoionization. Stasinska & Leitherer (1996), 
for example, restricted their analysis to metal-poor {Z < 
Zq/A) objects. As a matter of fact, statistical proper- 
ties of starbursts are not definitely established. A metal- 
dependent IMF was early proposed as an explanation for 
the [Oiii]a5oo7/H/3 decrease at high Z (Terlevich 1985; 
Shields & Tinsley 1976). Although a standard IMF seems 
in agreement with the bulk of observations (Leitherer 
1998), the IMF slope (see e.g. Eisenhauer et al. 1998; Greg- 
gio et al. 1998) and mass cut-offs (Goldader et al. 1997) 
are still under debate. 

Our aim is hereafter to study the evolution of the spec- 
tral properties of aging starbursts with the help of models 
taking into account the effects of mctallicity, geometry, 
dust and the excitation level of the gas. To avoid normal- 
ization problems, we selected relative properties indepen- 
dent of distance (colors, equivalent widths and line ratios). 
We focused on emission line ratios at close wavelengths - 
thus reducing reddening effects - to study the ionizing 
spectrum and the excitation of the gas. Equivalent widths 
and colors are used preferentially to study age effects, as 
well as the impact of the spatial distribution of stars, gas 
and dust. 

The observational sample is presented in Sect. 2. The 
coupling of the codes PEGASE, to model the evolution 
of star formation and stellar emission, and CLOUDY, for 
a consequent photoionization of the gas by massive stars 
in a given geometry, is described in Sect. 3. The relation 
between the mctallicity Z and the ionization parameter U 
is analyzed in Sect. 4, while equivalent widths and colors 



are considered in Sect. 5. The contribution of an under- 
lying population is analyzed in Sect. 6. Discussion and 
conclusion are respectively in Sect. 7 and 8. 

2. A selected sample of starbursts and Hll 

galaxies 

Our data set (Table 1) is a selection of 754 objects from ob- 
servational samples of galaxies having properties charac- 
teristic of starbursts: intense emission lines (French 1980; 
Veilleux & Osterbrock 1987; Terlevich et al. 1991), strong 
far-infrared emission (Leech et al. 1989; Veilleux et al. 
1995), or excess of blue or ultraviolet emission (Storchi- 
Bergmann et al. 1995); in addition, Contini et al. (1998) 
require galaxies to have a bar. Objects dominated by an 
active galactic nucleus (AGN) are excluded from our sam- 
ple, though some of them may be faintly contaminated by 
nuclear activity. Effects of aperture and differences in data 
reduction processes are considered in Sect. 2.2. 

2.1. Spectral classification of the original sample 

Objects listed in the source papers are classified in two 
families. The first one includes metal-poor Hii galaxies 
(H llGs) and extra-nuclear H ii regions, both likely devoid 
of active nucleus. The second one gathers galactic central 
regions where the ionized-gas emission could be partly due 
to shocks (Kim et al. 1998) or AGNs (see Heckman 1991 
and references therein). To select only starburst nucleus 
galaxies (SBNGs) in this last class, we keep the classi- 
fication of the authors, which is based on the pioneering 
works of Baldwin et al. (1981, hereafter BPT) and Veilleux 
& Osterbrock (1987, hereafter VO). 

The Terlevich et al. (1991) sample is mainly composed 
of HiiGs with some possible SBNGs. The authors used 
two of the BPT criteria based on the [Oiii]a5oo7/H/3 and 
[Oii]a3727/[Oiii]a5oo7 ratios. We also analyzed their sam- 
ple with the VO method and obtained the same classifi- 
cation. 

Contini et al. (1998) adopted the VO criteria. Their 
sample includes SBNGs and giant Hii regions located 
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Line 


PEGASE+CLOUDY 


Garcia- Vargas ct al. (1995) 


log(H/3) [ergs '-] 


38.83 


38.84 


[Oii]a3727/H/3 


3.10 


3.07 


[OIII]a5007/H/3 


0.31 


0.25 


[Oi]a6300/H/3 


0.04 


0.04 


[Nn]AG584/H/3 


1.34 


1.34 


[Sii]a(,71(,/H/? 


0.57 


0.58 


[Sll]Ai,7:u/n,:i 


o.;i9 


0.10 


[SIII]a9069/H/3 


0.36 


0.34 



Line 


PEGASE + CLOUDY 


Stasinska & Leitherer (1996) 


log(H/3) [ergs-^] 


40.56 


40.56 


[On]A3727/H/3 


2.57x10"' 


2.75x10"' 


[Oni]A5oo7/H/3 


8.18x10"' 


7.37x10"' 


[Oi]a6300/H/3 


2.60x10"^ 


3.13x10"^ 



Table 2. Comparison with predictions of coupled models available in the literature. 
Top: Garcia- Vargas ct al. (1995); age = 1 Myr, Z = Zq, log(J7) = -3.11. Bottom: 
Stasinska & Leitherer (1996); age = 1 Myr, Z = Zq, log(C/) = -2.60. 



in IRAS barred spiral galaxies, selected from Mazzarella 
& Balzano (1986) and listed in the Lyon-Meudon Extra- 
galactic Database (LED A). 

Two samples contain only SBNGs located inside IRAS 
galaxies. Veilleux et al. (1995) followed both BPT and 
VO criteria, while Leech et al. (1989) used only some of 
the criteria of BPT. Actually, an important part of their 
sample lies far beyond the limit reported by VO between 
Hll-like regions and AGNs in the [Oiii]/H/3 vs. [Oi]/Ha 
and [Oiii]/H/3 vs. [Nii]/Hq! diagrams. These objects are 
clearly inisclassified and have been excluded from our se- 
lection. Storchi-Bergmann et al. (1995) reported observa- 
tions of emission line galaxies of various types. Following 
Coziol ct al. (1998), all compact galaxies are hereafter clas- 
sified as H II galaxies, whether dwarf or not. Finally, we 
include in our selection the observations of five starbursts 
from the sample of Balzano (1983) reported by VO, and 
the data of French (1980) on 14 HllGs. 

2.2. Homogeneity of reddening, aperture and absorption 
corrections 

Most line ratios considered here are computed from close 
lines; so, reddening effects should be small and will be 
neglected in the following. 

In most samples, the contamination, through large 
apertures, of the starburst light by the host galaxy pop- 
ulation is likely to increase the continuum emission and 
to dilute the equivalent width of emission lines. To avoid 
this problem, we will restrict the analysis of colors and 
equivalent widths to the sample of Contini et al. (1998). 
These data were acquired with a long-slit spectrograph, 
but the spectra presented by the authors correspond to 
Ha emitting regions exclusively. Hence, the contamination 



of the starburst continuum by the environment of the host 
galaxy should be very weak. 

The host galaxy can also modify the apparent line ra- 
tios involving Balmer lines through the presence of stellar 
absorption lines. Emission line fluxes are corrected for this 
effect as follows: 

pcorr jTiobs / -i _i_ line 1 /i \ 

^line - ^line [^^ + J . (1) 

where i^i-"" and are respectively the absorption- 

corrected and the apparent emission line fluxes, and W^^l 
and are the absorption and measured emission 

equivalent widths. Not all the samples were corrected (Ta- 
ble 1). The correction should be negligible for the Ha line 
because Wg^" < Wg^". It is stronger for H/3 (up to 100% 
in Leech et al. 1989) and the ratio [Oiii]/H/3 may be af- 
fected. This ratio is systematically higher in the samples 
of French (1980) and Terlevich et al. (1991), mainly com- 
posed of HIIGs, than in the other ones, which include 
mainly SBNGs. Is this trend real, or is it due to the fact 
that these two data sets are not corrected for underlying 
stellar lines? Note that Storchi-Bergmann et al. (1995) al- 
ready noticed that [O iii]/H/3 is higher in H II galaxies than 
in SBNGs. Moreover, the correction for absorption lines 
should be very small in H llGs, since these objects have 
simuhaneously large VFgJ?/ (67 A in Terlevich et al. 1991) 
and low M^h/s'^. Hence, we conclude that the trend we ob- 
serve is real and not due to the discrepancies between the 
reduction procedures. 

3. The coupling of spectral evolution and 
photoionization 

The evolutionary synthesis code we use, PEGASE, takes 
into account metallicity and dust effects. PEGASE com- 



4 



E. Moy, B. Rocca-Volmerange, M. Fioc; Photoionization and star formation in starbursts 



Element X 


(X/H)z0 


(X/H)z 


He 


0.107 


0.08096 + 0.01833 x Z/{Q.7Zq) 


B 


2.63x10"^^ 


(B/H)zp, X Z/Zq 


C 


3.63x10 


(C/H)zp, X ZIZq 


O 


8.51x10'" 


(0/H)zp X ZIZq 


F 


3.63x10^'* 


(F/H);,^ X ZIZq 


Na 


2.14x10'" 


(Na/H)z„ X ZIZq 


P 


2.82x10"' 


(P/H)zq X ZIZq 


CI 


3.16x10'' 


(Cl/H)z0 X ZjZr^ 


Ar 


3.63x10'" 


(Ar/H)z„ X ZIZq 


Fc 


4.86x10'" 


(Fe/H)z<, X ZIZq 



log(N/0) 

12.43 - 3.77(log(0/H) + 12) + 0.26(log(O/H) + 12)^ 



Table 3. Adopted solar abundances and prescriptions for abundance 
variations with Z 



putes the stellar spectral energy distributions (SEDs) and 
the metallicities of starbursts and galaxies of the Hubble 
sequence at any stage of evolution, within the metallicity 
range Z = lO"* to 10' ^ Typical parameters of PEGASE 
are the star formation rate (SFR) and the initial mass 
function. The new version]^ used hereafter, PEGASE. 2, is 
based on the evolutionary tracks of Girardi et al. (1996), 
Fagotto et al. (1994a,b,c) and Bressan et al. (1993). The 
AGB to post- AGE phases are computed following the pre- 
scriptions of Groenewegen & de Jong (1993) models. The 
synthetic stellar spectral library is from Kurucz (1992), 
modified by Lejeune et al. (1997) to fit observed colors. 
For details, see Fioc & Rocca-Volmerange (1999; 2000, in 
preparation) 

The photoionization code CLOUDY (version 90.04, 
Ferland 1996) predicts the spectra of low- to high-density 
astrophysical plasmas in the most extreme astrophysi- 
cal sites, such as starbursts and quasar environments. It 
takes into account recent changes in atomic databases and 
new numerical methods. In addition to its physical per- 
formances, the code was chosen for its easy Web access, 
its documentation and its friendly on-line help. CLOUDY 
proposes options to explore the influence of geometrical 
factors: type of geometry - spherical or plane parallel -, 
covering factor ri/(47r), filling factor /, and the distance 
R from the source to the plasma. 

The coupling of PEGASE. 2 with CLOUDY allows to 
compute coherently the stellar and nebular emission, and 
to fit simultaneously the continuum and lines from stel- 
lar populations embedded in a gas cloud. The metallicity 
Z'(i), and the stellar SED - in particular the energy distri- 
bution of the ionizing photons - provided by PEGASE are 

^ This code is available at 

http://www.iap.fr/users/fioc/PEGASE.html or by anony- 
mous ftp at ftp.iap.fr in /pub/from_users/pegase. 



used as inputs by CLOUDY. Hereafter, we call Qh" the 
emission rate of ionizing photons. This quantity is used to 
normalize the SEDs computed by PEGASE. 

The results of the current version of the code have been 
compared to those obtained by similar "coupled models" . 
They are presented in Table 2 for two sets of predictions: 
the model 16 from Garcia- Vargas et al. (1995b), and the 
^Q-model computed by Stasihska & Leitherer (1996) for 
a 10^ Mo-cluster. In both cases, we used inputs matching 
the reference model (IMF, age, geometrical parameters 
and chemical abundances). For the first comparison test, 
the agreement is very good (Table 2, top). Our results also 
agree relatively well with those of Stasinska & Leitherer 
(1996), given that the evolutionary tracks, the spectral 
libraries and the photoionization codes are different in the 
two models. 



3.1. Star formation parameters 

An instantaneous star formation episode ("pure star- 
burst") is assumed for starbursts - a scenario well suited 
to study emission line properties. In Sect. 6, we will 
also consider an underlying stellar population formed ei- 
ther during a previous burst or continuously. We tenta- 
tively tested several IMFs of the form dN/dm oc to~", 
m e [Afiow, A^up]: a = -2.35, m e [0.1, 120] (Salpeter 
1955); a = -3.35, m G [0.1, 120] Mq; and a = -2.35, 
m e [3,30]Mq. The last one is truncated, as suggested 
by various infrared studies (Rieke et al. 1993; Langon & 
Rocca-Volmerange 1996). The initial metallicity is left as 
a free parameter. 
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Fig. 1. Models of pure instantaneous starbursts for various filling factors (/) and metallicities (Z). 
The data are from Terlevich et al. ("plus" signs), Veilleux & Osterbrock (asterisks), French (filled 
circles), Veilleux et al. (diamonds), Contini et al. (triangles), Storchi-Bergmann et al. (squares) 
and Leech et al. (crosses). 



3.2. Physical properties of the gas 

The emission line spectrum of an ionized nebula depends 
on the combination of the ionizing spectrum, the chem- 
ical composition of the gas and the so-called ionization 
parameter U. 

The elemental abundances of the gas have been made 
consistent with the metallicity of the ionizing stars. The 
solar abundances used here are given in Table 3. All the 
abundances were scaled linearly according to the metal- 
licity, except for nitrogen and helium. For the likely sec- 
ondary element N, we adopted the law proposed by Co- 
ziol et al. (1999), multiplied by a factor 1.5 to account 
for the excess of nitrogen observed in starbursts. For He, 
we followed a prescription proposed by Dopita & Kewley 
(private communication; see Table 3 for details). 



The ionization parameter U is defined as: 



U 



47ri?^riHc' 



(2) 



where nu is the hydrogen density, c the speed of light, and 
Rs the Stromgren radius (i.e. the radius of the ionization 
front). Taking into account the inner radius R (e.g. Evans 
& Dopita, 1985), Rg is roughly equal to: 



R. 



3Qh" 



R' 



1/3 



(3) 



where cjb is the case B temperature-dependent recom- 
bination coefficient and / is the filling factor. This for- 
mula highlights the impact of each parameter on U, but is 
only approximate, since it assumes a constant temperature 

through all the nebula. The effective Stromgren radius is 
obtained only at the end of a CLOUDY calculation. 

The spherical geometry (Stasihska & Leitherer 1996; 
Gonzalez-Delgado et al. 1999) and plane-parallel geom- 
etry (Garcia- Vargas & Diaz 1994; Garcia- Vargas et al. 
1997) are extreme situations for which the right-hand 
side of eq. (3) is respectively dominated by the first and 
the second term. The adopted geometry was spherical. 
The hydrogcm density nu was 100 cm~^ in most mod- 
els, but we also computed a few cases for nn = 10 and 
riH = 1000 cm~^, to assess the specific impact of den- 
sity on emission line ratios. We set the number of ionizing 
photons to Qua ~ 10^^ photons s~^, corresponding to a 
maximum luminosity of ~ 10^*^ ergs""'^ in the Ha emission 
line - a value characteristic of intense star formation sites 
(Shields 1990). For emission fines, the exact value of Quo 
is not important as long as the density does not exceed 
the critical limit of collisional de-excitation, below which 
models with the same U, ionizing spectrum and metal- 
licity predict the same line ratios. To obtain reasonable 
values for U, we varied the volume filling factor between 
10~^ and 1 and set Rto 4 pc. 
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Fig. 2. Same as Fig. |l] for 

dening. 

We assume that the gas is distributed in a spherical 
shell covering a solid angle fl around the star cluster. The 
covering factor 17/ (47r) was allowed to vary in the interval 
[0.1, 1]. The continuum luminosity from the starburst F\ 
(ergs~-^ A~^) at any wavelength A is: 

Fx = ^"'^ + F^"^ + Ff °". (4) 

The transmitted stellar continuum F^''^'^^ and the nebular 
continuum F^°^ both scale with r2/(47r). F^^'^^^, scaling as 
(1 — n/{4:Tr)), refers to the part of the stellar spectrum 
not crossing the gas (and thus not contributing to the 
ionization process). 

3. 3. Extinction of the spectrum 

Assuming a dust screen at the outer edge of the H ii region, 
the continuum flux is reddened as follows: 

Fx (Ff '^"^^ + F^'='^)e-^J'''' + Ff ^"e-^S*°" . (5) 

Emission line fluxes are also multiplied by e""^^ . The 
optical depth Ty^"" was a free parameter, as well as the 
^stciy^ncb j.g^^jQ (■< \y This option takes into account a 
possible dust concentration outside the nebula rather than 
on the line-of-sight not covered by the gas. The variation 
of the optical depth with wavelength is derived from the 
metallicity-dependent extinction laws of the Magellanic 
Clouds and our Galaxy. 

4. Results on emission line ratios for 
instantaneous starburst models 

We have explored the impact of the IMF, the age of the 
starburst, the metallicity and the ionization parameter on 
the following optical line ratios: [O iii]a4363/[0 iii]a5007j 
[Oi]a630o/Hq;, [Sii]AA67i6,673i/Ha, [N ii] A6584/Ha and 
[O iii]a5007/H/3. The value of / was varied between 10^^ 
and 1. The corresponding log(C/) (as defined in Sect. 3.2) 
belong to [—5.5,-1.5]. Note that similar results could also 
be derived by varying the inner radius R (in a plane- 
parallel geometry), the density tih or a combination of 
both. 




-2 - ■ - 

Log ( [SIN] X 9069 / [Sll] AX 6716,6731 ) 



line ratios sensitive to red- 



Metallicity effects 

Emission line ratios are sensitive to the metallicity 
through two competing effects. On the one hand, the opac- 
ity increases with Z, lowering the stellar effective temper- 
atures and thus the hardness of the ionizing spectrum. A 
trend towards lower excitation levels is therefore expected 
in high-Z starbursts, even if there is no dependence of the 
IMF or the geometry on the metallicity. This effect is re- 
inforced at high Z by the efficient cooling by metals and 
thus the low electronic temperature To- 

On the other hand, the energy in metal emission 
lines tends naturally to increase with higher elemental 
abundances. For oxygen lines, this process dominates at 
Z < 1/5^0, while the opposite happens at higher metal- 
licity (e.g. McGaugh 1991). As a consequence, two models 
with metallicities respectively below and above this limit 
can predict comparable emission lines. For this reason, 
we present hereafter our predictions for Z in the range 
[1/50, 3/2]Z0 if no degeneracy happens in the considered 
diagram, but in [1/5,3/2]Zq otherwise. 

4-.2. Effects of the geometrical filling factor f 

For a given density njj, varying the filling factor / changes 
the optical depth per unit length and the ionization struc- 
ture: the smaller /, the smaller the mean ionization level, 
and the lower the luminosities of high excitation lines such 
as [Oiii]a5007- Low-excitation lines behave inversely. 

Models of pure instantaneous starbursts, with metallici- 
ties in the range [l/5,3/2]Z0, are compared to the data 
in Fig. |l| for six values of the filling factor /: 10~^, lO"'', 
IQ-'^, 10"^ 10-\ and 1. The slope of the IMF is -1.35, for 
a mass range of 0.1-120 M©. The effects of dust are not 
taken into account in the models plotted on Fig. g. The 
data are largely covered by the models when at least two 
parameters vary - here chemical and geometrical. Inter- 
estingly, different types of objects seem to occupy different 
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Fig. 3. Comparison of the data with five sequences obtained by varying the relative weights of the 
high (/ — 1) and low (/ = 10"^) excitation components. Solid diamonds are the points for which 
the low excitation component contributes to 0, 1, 10, 50, 90, 99 and 100% of the H/3 flux (from 
right to left). 
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Table 4. Variations of line ratios induced by a change of density. The 
age is OMyr and the metallicity is Z^jh. 



zones in the Z-\J plane: Hii galaxies are in better agree- 
ment with low-Z (1/50 Zq to 3/4 Zq), high-[/ models 
(/ G [10-'^, 1]); SBNGs are compatible with low values of 
/ (10"^ to 10-^), but encompass a wide range of metal- 
licities. 

We also report on Fig. |^ the comparison of predictions 
with observations in following diagrams: [Oii]a3727/H/3 
vs. [Oiii]a5007/H/3 and [Oii]a3727/[Oiii]a5007 vs. 
[Siii]a9069/[Sii]aa6716,673i- Notc that the [Oii]/H/3, 
[Oii]/[Oiii] and [Siii]/[Sii] ratios are sensitive to redden- 
ing and, for this reason, we do not include them in the 
discussion. We can only notice that, given the uncertain 
impact of dust on these ratios, our results are compatible 
with the observations. 

Jf..3. The apparent relation between Z and U: degeneracy 
problems 

A mean U-Z relation fitting the data could easily be de- 
rived from Fig. 0, as in Dopita & Evans (1986). However, 
this relation would suffer from too many uncertainties: 
first, it would not take into account the high Z-low Z 



degeneracy; second, age effects are very similar to those 
induced by the metallicity. This is illustrated in Fig. |[ 
where [Oi]/Ha and [Sii]/Ha vs. [O iii]a5007/H/3 are plot- 
ted for Zq/5 at < = 0, 4, 5 and 6Myr. Variations of the 
IMF give rise to similar uncertainties, as shown on Fig. |[ 
the lower number of massive stars for a steeper or trun- 
cated IMF causes a strong decrease of both [O ill] /H/3 and 
[Ol]/Ha, just as if the metallicity was higher. 

Due to this age-metallicity-IMF degeneracy, a mean- 
ingful Z-U analytical relation cannot be established on 
the basis of line ratio diagrams. Moreover, Fig. |l| clearly 
shows that the dispersion of U, rather than the mean value 
of this parameter, is linked to the metallicity: while star- 
bursts exist at any U at low Z, only \ow-U starbursts are 
observed at high metallicity. We consider a possible cause 
of this trend in the next paragraph. 

4-.4- The upper envelope of the data sequences 

In the [Oi]/Ha vs. [O ill] /H/3 and [Six] /Ha vs. 
[Oiii]/H/3 diagrams, many data points lie above predic- 
tions, even when the misclassified objects of Leech et al. 
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Fig. 5. Time evolution of the [O iii]a5oo7/H/3 ratio for 
instantaneous starbursts. The decrease of the ratio 
after 5 Myr is very rapid and puts severe constraints 
on the age of starbursts (see text for details) 



(1989) are excluded from the sample. To check whether 
this discrepancy can be solved by varying the density, we 
have performed two series of calculations at Z ^ Zq/5 
with TT-H = 10 and 1000 cm~^. We chose this metallicity 
because the corresponding sequence defines the upper en- 
velope of our model grid in these diagrams. The results are 
shown in Table ^. The main one is that model sequences 
with different densities largely overlap: clearly, the outliers 
cannot be explained by a variation of nn- 

Many authors have proposed that a faint contribu- 
tion from shocks or AGNs could explain the intensity 
of the [Oi] line in starbursts (see e.g. Stasiriska & Lei- 
therer 1996). We consider here an alternative possibil- 
ity where there are, inside the starburst, both high- and 
low-excitation zones, due for example to a distribution of 
gas and stars more complex than the classical geometry 
adopted in our calculations. To test this hypothesis, we 
computed a series of sequences obtained by varying the 
relative weight of high- and low-excitation models (/ = 1 
and / = 10~^, respectively) for five metallicities. The re- 



sults (Fig. ^) are in much better agreement with the data 
than pure U sequences, as they explain both the trend 
and the dispersion of the observations. Another advantage 
of this explanation is that no other source of ionization, 
such as AGNs or shocks, is required. If this hypothesis is 
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Fig. 6. Impact of the IMF on [Oi]/Ha vs. [Oiii]/H/3 
for two metallicities {Zq/5 and Zq). In each case, 
three / sequences are plotted: A (a = —2.35, Miow = 
0.1 Mq and Mup = 120 M©), B (a = -3.35, Afiow = 
0.1 Mq and Mup = 120 M©) and C (a = -2.35, 
Miow = 3 Mq and M^p = 30 Mq). 



correct, the evolution of the apparent U with Z in fact 
corresponds to the presence in low-metallicity objects of a 
high-J7 component, which is absent in high-Z starbursts. 
We come back to this point in the discussion. 

4-5. Constraints on age and IMF 

Fig. H shows the evolution, up to 6 Myr, of the 
[O iii]a5007/H/3 ratio for the models plotted on Fig. |l| and 
^. The constraints on age shown by this plot are very 
strong, though noticeably different for Hii galaxies and 
SBNGs: Hii galaxies, of high-excitation level, need the 
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Fig. 7. Pure starburst models for various ages (0 to 6Myr) and Ty between to 6. The ratio r^^'^'V^r'' = 0-5 
is equal to 1. The covering factor is 0.1. Data (triangles) are from Contini et al. (1998). The solid lines indicate 
the predictions of the 0/(47r) = 0.1 model grid. The two dashed lines show the t = OMyr (upper hne) and 
t — ^ Myr (lower line) model results for / (47r) = 1 . 



high energy photons emitted by massive stars, so that only 
young (<3Myr) models are compatible with observations 
(see Sect. 7); SBNGs, requiring a lower excitation level, 
span a wider range of ages. In any case, the rapid drop 
of the [O III] /H/3 ratio after 5 Myr requires the presence of 
stars younger than 6 Myr. 

Line ratios depend on the hardness of the spectrum, and so 
on the IMF. The [O i] line, emitted from the partly ionized 
zone, requires high-energy photons and is thus particularly 
sensitive to the high-mass end of the IMF. Fig. ^ presents 
the sensitivity of [O i] a6300 /Ha and [O ill] /H/3 to the initial 
mass function. A Salpeter IMF is favored by our results, 
in good agreement with recent studies (e.g. Stasihska and 
Leitherer 1996; Garcia- Vargas et al. 1995). An interesting 
point is that the steeper and truncated IMFs, while clearly 
excluded by the data for high-[/ HiiGs, are acceptable for 
metal-rich starbursts. 

5. Results on equivalent widths and colors: 
covering factor and extinction 

In the following, we aim to keep our previous fits of emis- 
sion line ratios and, simultaneously, to reproduce M^(Ha) 
and B — V and V — R colors from the sample of Contini 
et al. (1998). The models plotted on Fig. [l] and | pre- 
dict Ha equivalent widths above 500 A, while the observed 
tl^(Ha) seldom exceed 300 A. Predicted equivalent widths 
in excess to the observations are a well-known problem 
(Bresolin et al. 1999 and references therein). The hypothe- 
sis of a "density bounded" nebula (Leitherer et al. 1996) is 
ruled out by the observed intensities of the low-ionization 
lines; the reduction of the thickness of the partly-ionized 
zone would lead to lower [Sll]/Ha and [Oi]/Hq! than in 
Fig. 0, and the models would not fit any more these line 
ratios. 



A low covering factor ri/(47r) reduces the equivalent 
widths by increasing the stellar emission relative to the 
nebular one, without changing line ratios. A similar expla- 
nation of the faint W{llf3) of blue compact and irregular 
galaxies was suggested by Mas- Hesse & Kunth (1999). 

Fig. 1^ shows the influence of dust and age in the 
WiRa) vs. B ~V and W{B.a) vs. V ~ R planes, for 
r2/(47r) = 0.1. For comparison, we also show t — and 
t — 6 Myr sequences for il/(47r) — 1. We plot only Z = Zq 
models, since the metallicity has only weak effects on 
B — V, V — R and VF(IIq!). Similarly, for such a low value 
of ri/(47r), / hardly affects the colors and the only mod- 
els on Fig. are for / = 1. Ages between and 6 Myr, 
^stoiy^gas _ ^ ^gas range [0,6] are suggested 

by the VF(Ha) vs. B — V diagram. However, ry > 4 is 
excluded by the V ~ R data, and ages > 6 Myr are incom- 
patible with line ratios (see also Sect. 7.2). 

A ratio rf •^^Vrf''" ~ 0.5 (following Calzetti et al. 1997; 
see also Fanelli et al. 1988; Keel 1993; Mas-Hesse & Kunth 
1999) can simultaneously fit B — V and V — R, but requires 
very high extinctions, up to Ty"^ = 8 (Fig. while Con- 
tini et al. (1998) derived a maximum value of 4 from 
the Balmer decrement. The B — V oi the reddest objects 
are considered as anomalous (Contini, private communi- 
cation) and are not reproduced by any model. 

6. The contribution of an evolved population 

The distribution of the data in the VF(Ha) vs. B — V or 
V — R diagrams is very difficult to explain with pure in- 
stantaneous starburst models. The main problem is the 
very low covering factor required. A value of 0.1 seems 
in contradiction with previous estimates of the escaping 
fraction of ionizing ohotons in HII regions (see Sect. 7.2). 
An alternative hypothesis, able to explain the equivalent 
widths and colors of the Contini et al. (1998) sample with 
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Table 5. Parameter values of starburst + host scenarios, (a): Host galaxy age (in Myr) 
when starburst occurs, (b): star formation parameter v (in Myr~^) of the host galaxy, (c) 
Time-scale of gas infall onto the host galaxy (Myr) . (d) : fraction of ionizing photons emitted 
by the host population when the starburst starts in models A, B, C, D and E (see text). 



the much more reasonable value n/(47r) = 0.5, is the pres- 
ence of evolved stars inside the starburst. Such a popula- 
tion, created during past formation episodes, can account 
for the observed weakness of equivalent widths by increas- 
ing the continuum (e.g. McCall et al. 1985; Diaz et al. 
1991). It may also explain the red colors of the sample 
without requiring huge amounts of dust. 

We considered two possible underlying populations. 
The first one consists of an instantaneous burst which oc- 
cured 100 Myr ago. This old burst is assumed to have been 
ten times more intense than the current one. The predic- 
tions of these models are plotted on Fig. ^. As in Fig. ||, 
T|*'^"/rj^^'^^ is set to 0.5. However, these models face two 
problems: first, it is still necessary to extend the age of 
the young burst to i = 6 Myr, unless the covering fac- 
tor is much lower than the value adopted on Fig. ^ (0.5); 
second, the sum of the evolved and the young population 
does not reproduce the colors. To reconcile this type of 
models with the data, it would be necessary to assume 
that the old burst was systematically ^ 100 times more 
intense than the current one. 

The second type of underlying population is assumed 
to have formed continuously. We computed a series of sce- 
narios matching the observations of standard Hubble se- 



quence galaxies at z = 0, and added a starburst. For the 
sake of consistency, the initial metallicity of the burst is 
the current metallicity Z(t) of the host galaxy and is pro- 
vided by PEGASE. 

The input parameters are chosen so that the starburst 
metalhcities belong to the range [1/5Zq, 3/2^©] derived 
from line ratios. Each scenario is defined by three param- 
eters: the age tgai of the host galaxy when the starburst 
occurs; the infall timescale tc used to parameterize the 
accretion rate of the galaxy; the gas-to-star conversion ef- 
ficiency, ly, relating the star formation rate ({t) to the gas 
density a{t) by the Schmidt law C{t) = va{t). 

Five models (A, B, C, D, E) were computed for each 
scenario, corresponding to five different contributions of 
the evolved population to the total emission: A is the pure 
host galaxy, E is the pure starburst and B, C and D are 
intermediate cases. Note that the star formation does not 
stop after the burst, but returns to the quiescent regime. 
The fraction of ionizing photons due to the underlying 
population immediately after the burst is given in Table 
5, as well as the values of tgai, ic, and v for the considered 
scenario. 



Fig. 10 presents the impact of the parameters on colors 
and equivalent widths. We compare the results of model 
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D at Z = Z0/2 to the sample of Contini at al. (1998). 
We also show the results of model C for Z = 1/2Zq, and 
the predictions of model D for 3/2^0. The filling factors 
are lO'^ at 1/2^0 and 5 x IQ-^ at 3/2^0. As in Fig. |, 
^steiy^neb jg g_g^ Howcver, the covering factor (0.5) 

is significantly higher, and the Ty°^ range (0 to 4) is much 
more reasonable. 

Another important point is that, contrary to the pure star- 
burst case, ages of 6 Myr and more are compatible with 
emission line ratios. This is due to the presence of massive 
stars, formed after the burst in the underlying population. 
The presence of such stars ensures that the line ratios are 
virtually unchanged compared to their values during the 
burst. Contrary to the case of a burst that occured 100 
Myr ago, these scenarios are able to explain simultane- 
ously the colors, the equivalent widths and the line ratios 
in our sample. They do not require very high levels of 
extinction, and, maybe more important, do not imply a 
covering factor lower than 0.5. For all these reasons, the 
hypothesis of an underlying population formed continu- 
ously is our preferred one. 

7. Discussion 

7.1. A physical link between the ionization parameter and 
the metallicity? 

Dopita & Evans (1986) proposed an analytical relation be- 
tween the metallicity and the ionization parameter to de- 
scribe the observational distribution of line ratios in their 
sample of extragalactic Hii regions. In Sect. 4.3, we saw 
that such a relation suffers from many degeneracies due 
to age, IMF and metallicity. Moreover, it seems that the 
quantities which are really related are the dispersion of U 
and the metallicity. 

What could be the physical link between Z and Ul Ac- 
cording to one hypothesis, this is a link through geomet- 
rical effects. Following Castor et al. (1975) and assuming 
their equations are valid also for a star cluster, the inner 
radius i? of a gaseous shell surrounding a massive star or 
cluster is: 

i?oc(^:^^ , (6) 

where E is the energy deposition rate of a star (or cluster) 
of age t in the interstellar medium (ISM) of density p. 
Leitherer et al. (1992) found that E is constant for an 
instantaneous burst until ~ 6 Myr and scales linearly with 
Z, regardless of the age. Hence, high-Z H 11 regions expand 
faster. 

We therefore propose the following scenario to explain 
the distribution of starbursts in the line ratio diagrams. 
For some reason, a site of intense star formation appears 
in a galaxy with a metal-poor interstellar medium. This 
starburst can be composed of many individual H 11 regions 



surrounding star clusters. As these sub-components evolve 
with age, the Hii regions expand because of stellar winds. 
At the same time, the enrichment of the ISM by massive- 
star ejecta begins. If star formation goes on, very young 
components will coexist inside the starburst with more 
evolved components. The former will dominate the nebu- 
lar emission during the first few Myr, and the latter after. 

Two additional effects can explain the lack of high-.Z 
and high-C/ objects: first, at high metallicity, individual 
Hii regions expand faster, leading to a lower ionization 
parameter; second, according to this scenario, the bulk 
of nebular emission in high-Z starbursts is produced by 
evolved, low-excitation H 11 regions. Checking the validity 
of this hypothesis requires to model also the expansion of 
nebulae. This is the subject of a forthcoming paper. 

7.2. The equivalent-width problem 

The problems encountered by pure starburst scenarios to 
account for the weakness of the Ha equivalent widths from 
the sample of Contini et al. (1998) are several. The main 
one is that a very low covering factor is required. From 
Fig. 1^ and ||, we deduce that a value r2/(47r) ~ 0.1 is nec- 
essary to reproduce both l/F(Ha) and the line ratios in the 
case of a pure starburst. This means that 90 % of the Ly- 
man continuum photons do not contribute to ionization, 
but are directly absorbed by dust or escape from the star- 
burst. This is in contradiction with many previous studies 
using Ha or H/? luminosities to estimate the mass of stel- 
lar clusters inside starbursts. In many cases, the agreement 
between such estimations and those based on ultraviolet 
continuum luminosities is excellent (Gonzalez-Delgado et 
al. 1999). Some discrepancies appear frequently (Leitherer 
et al. 1996; Vacca et al. 1995), but they are however weak 
enough to exclude an absorption fraction of the ionizing 
photons by the gas as low as 10%. Moreover, observa- 
tions of diffuse Ha emission in nearby galaxies lead to an 
estimation of 15-50 % for the fraction of ionizing photons 
escaping from individual H 11 regions (Kcnnicutt 1998 and 
references therein). 

Such a systematic low value of the covering factor is 
therefore clearly ruled out by the observations. Note also 
that pure starbursts models, even with 17/(47r) = 0.1, can- 
not reproduce simultaneously the emission line ratios and 
the equivalent widths (see Sect. 5). Hence, the presence 
of an evolved stellar population coexisting with a star- 
burst, as analysed in Sect. 6, seems mandatory to explain 
equivalent widths and colors, even for the data of Con- 
tini et al. (1998), who strictly limitated their apertures to 
Ha-emitting areas. A similar conclusion for galaxies ob- 
served through 5-arcsec. apertures was drawn by Langon 
& Rocca-Volmerange (1996) from the near-infrared spec- 
tral synthesis of starbursts. 
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Fig. 9. Ha equivalent widths and colors for instantaneous + lOOMyr-old starbursts. The 
metallicity of the young starburst is solar. The old burst is assumed to have been ten times 
more intense than the current one. The model sequences are labeled in terms of the age of 
the present starburst and of the value of Ty°'^. 




Fig. 10. Ha equivalent widths and colors for instantaneous burst + continuous star for- 
mation scenarios. The thick solid lines show the predictions of model D for different ages 
and extinction levels of the burst. The metallicity of the burst is 1/2 Zq. The dashed lines 
indicate the shift of the right part of the grid when the contribution of the underlying pop- 
ulation is increased (model C). The thin solid lines show the shift of the low part of the grid 
when the metallicity is equal to 3/2 Zq. 



7. 3. Nature of the underlying population 

Including underlying populations in the models allows to 
fit the data in a satisfactory way in the PF(Ha) ms. B — V 
and V — R diagrams with ri/(47r) — 0.5, which means that 
~ 50% of the ionizing photons escape from the gaseous 
nebula. This value is close to the 60 % fraction suggested 
by Mas-Hesse & Kunth (1999). The nature of this popu- 
lation remains however elusive. 

We have tested an old (100 Myr old) instantaneous 
burst as a possible underlying population. The agreement 
with the data, at first sight, is relatively good (Fig. In 
this case, however, the age of the youngest burst has to 
extend up to 8 Myr to explain the faint end of the W(Ha) 
observational distribution. At this age, our models are in- 
compatible with line ratio data. Note that there is no Wolf- 



Rayet(W-R) star SED in our spectral library; instead, the 
spectra of the hottest stars available in the library are used 
during W-R stages. The implementation of W-R SEDs in 
the library could in principle extend the range of ages 
compatible with line ratios up to 6 Myr or more. As an 
example, Gonzalez-Delgado et al. (1998) obtained an age 
between 6 and 9 Myr for IRAS 08334-6517. However, this 
conclusion is based on models with an IMF truncated at 
= 30 Mq. Such a low value is clearly excluded for 
the bulk of our sample (Fig. ^). Moreover, the role of W- 
R stars in starburst nebular emission is still under debate. 
According to the most recent studies, the bulk of nebular 
emission is due to stellar populations younger than 3 Myr 
(Bresolin et al. 1999). 

We have assumed that the old burst was only ten times 
more intense than the current one. Changing this value to 
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one hundred, for example, would move the predictions on 
Fig. ^toward the bottom- right and might explain the bulk 
of the data with ages lower than, or equal to, 4 Myr. In this 
case, however, such a difference between past and present 
star formation remains to be explained. 

In this context, the models plotted on Fig. [l^ appear 
as the most attractive ones. Scenarios with an underlying 
population providing from 2.5 up to 20% of the number 
of ionizing photons are in acceptable agreement with the 
data and explain naturally all the observables (emission 
lines, colors and equivalent widths). The range of Ty is 
consistent with the Balmer decrement measurements of 
Contini et al. (1998). Moreover, the age and metallicity 
of the underlying population are typically those of normal 
spirals, suggesting that starbursts are normal events in 
galaxy evolution. These scenarios draw a coherent picture 
of the starburst phenomenon, and for this reason, remain 
the most plausible ones. 



8. Conclusion 

We have analyzed the emission line and continuum proper- 
ties of a large sample of starbursts and put new constraints 
on their stellar populations, metallicity, photoionization 
state and dust content. Our tool was the spectral evolu- 
tion code PEGASE coupled to the photoionization code 
CLOUDY. 

From the comparison of a dataset of emission lines 
with pure instantaneous starburst scenarios, we have been 
able to constrain the evolution of the ionization parameter 
U with the metallicity Z. At high abundance, the high- 
excitation component clearly detected at low Z is absent 
or very weak. A preliminary analysis based on energy de- 
position rate considerations is unsufficient to explain this 
trend, and a more refined analysis of the dependency of 
the H II region expansion with Z is needed. 

The second important result concerns underlying pop- 
ulations inside starbursts; they are detected even through 
small-aperture observations. This contribution is needed 
to reproduce simultaneously all the observables of our 
sample (line ratios, colors and equivalent widths). The 
underlying populations are simulated through the whole 
computation of the past star formation history of the 
galaxy hosting the starburst. In the scenarios favored by 
our analysis, the burst itself is supposed to be a brief, 
intense episode, preceded and followed by a quiescent 
regime. After the burst, this quiescent formation rate is 
sufficient to keep the line ratios compatible with the data. 

Future prospects should consider the problem of metal 
depletion on dust grains, since depletion changes the com- 
position of nebular gas. However, the question of the pres- 
ence of dust inside Hii regions will be solved only with 
more constraining data, e.g. those of the ISO satellite. Fi- 
nally, the large grids computed with the help of a coupled 
photoionization and stellar evolutionary synthesis model 



will be also used to study other samples with better statis- 
tics. 
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